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ABSTRACT 

In this paper we discuss the optical and X-ray spectral properties of the sources detected in 
a single 200 ks Chandra pointing in the Groth-Westphal Strip region. A wealth of optical 
Q-f photometric and spectroscopic data are available in this field providing optical identifica- 

q tions and redshift determinations for the X-ray population. The optical photometry and spec- 

$_i ■ troscopy used here are primarily from the DEEP2 survey with additional redshifts obtained 

^ \ from the literature. These are complemented with the deeper (r « 26 mag) multi-waveband 

C3 ■ data (ugriz) from the Canada France Hawaii Legacy Survey to estimate photometric redshifts 

and to optically identify sources fainter than the DEEP2 magnitude limit (Rab ~ 24.5 mag). 
We focus our study on the 2-10keV selected sample comprising 97 sources to the limit 
■ w 8 x 10 _16 ergs _1 cm -2 , this being the most complete in terms of optical identification 

5— i ' rate (86%) and redshift determination fraction (63%; both spectroscopic and photometric). 

We first construct the redshift distribution of the sample which shows a peak at z sa 1. This 
is in broad agreement with models where less luminous AGNs evolve out to z ss 1 with pow- 
erful QSOs peaking at higher redshift, z rs 2. Evolution similar to that of broad-line QSOs 
applied to the entire AGN population (both type-I and II) does not fit the data. We also ex- 
plore the observed Nh distribution of the sample and estimate a fraction of obscured AGN 
(Nh > 10 22 cm -2 ) of 48 ± 9 per cent. This is found to be consistent with both a luminosity 
dependent intrinsic Nh distribution, where less luminous systems comprise a higher fraction 
of type-II AGNs and models with a fixed ratio 2:1 between type-I and II AGNs. We further 
compare our results with those obtained in deeper and shallower surveys. We argue that a 
luminosity dependent parametrisation of the intrinsic Nh distribution is required to account 
for the fraction of obscured AGN observed in different samples over a wide range of fluxes. 

Key words: Surveys - galaxies: active - X-rays: galaxies - X-rays: diffuse background - 
cosmology: observations 



1 INTRODUCTION 

In the last few years the study of the diffuse X-ray background 
(XRB) has witnessed significant observational progress allowing 
detailed comparison with model predictions. The ultra-deep Chan- 
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dm surveys in particular, have demonstrated that most of the XRB, 
at both soft and hard energies, is resolved into discrete point sources 
(Brandt et al. 2001; Giaconni et al. 2002; Alexander et al. 2003), the 
vast majority of which are without doubt AGNs. To the first approx- 
imation, this finding has been a huge success for models that repro- 
duce the spectral properties of the XRB under the zero order as- 
sumption that it originates in a combination of obscured and unob- 
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scured AGN (Comastri et al. 1995; Gilli, Salvati & Hasinger 2001). 
Under more careful examination however, a number of inconsisten- 
cies emerge. Firstly, luminous (Lx > 10 44 ergs _1 ) heavily ob- 
scured type-II QSOs at z « 1.5 — 2, predicted in large numbers by 
the models, are scarce in the surveys above. Secondly, the redshift 
peak of the X-ray population lies below z = 1 in stark contrast 
with the model expectation of z « 1.5 — 2. 

Although these inconsistencies suggest that some revision of 
the models is almost certainly required (Hasinger 2003), observa- 
tional biases may complicate any interpretation. For example, about 
?s 25 per cent of the sources in the Chandra Deep Fields (CDF) 
are optically faint, R > 24 mag, rendering optical spectroscopy 
difficult or even impossible with current technology (Rosati et al. 
2002; Barger et al. 2003). Any information about the nature of these 
sources is therefore limited and they are proposed as best candi- 
dates for heavily obscured AGN (Alexander et al. 2001; Treister 
et al. 2004), likely to comprise a fraction the elusive population 
of high-z type-II QSOs. Moreover, the small field-of-view of the 
CDFs (0.07 deg 2 each) makes them sensitive to cosmic variance 
further complicating interpretation of the derived redshift distribu- 
tion. 

Wide-area shallower (« 10~ 14 erg s _1 cm~ 2 ) surveys are 
less affected by the observational biases above (e.g. Baldi et al. 
2002; Kim et al. 2004; Georgantopoulos et al. 2004). These sam- 
ples although of key importance, comprise a large fraction of unob- 
scured AGNs that are not representative of the sources responsible 
for the spectral shape of the XRB (r = 1.4; e.g. Gruber et al. 
1999). 

The evidence above suggests that deep surveys with relatively 
wide field-of-view are essential to improve our understanding of the 
XRB. Observational programs in this direction are already well un- 
derway such as the XMM-Newton Cosmic Evolution Survey [COS- 
MOS; 2 deg 2 , /x(0.5 - 2keV) « 5 x 10~ 16 ergs^cm- 2 ] and 
the Extended Chandra Deep Field South [E-CDF-S; Lehmer et al. 
2005; Virani, Treister & Urry 2006; 0.3 deg 2 , fx (0.5 - 2 keV) « 
1.1 x 10~ 16 ergs -1 cm -2 ]. In this paper we present results on a 
single 200 ks Chandra pointing, which is part of an on-going X-ray 
survey in the Extended Groth Strip region, which will eventually 
cover a total area of about 0.5 deg 2 to the depth above (200 ks per 
pointing). This sample, when completed, will be intermediate in 
terms of area coverage and depth to the CDFs and shallower wide- 
area surveys, minimising any observational biases affecting the 
ultra-deep fields and comprising a large fraction of obscured AGNs 
responsible of the XRB properties (Nandra et al. 2005). More- 
over, the Extended Groth Strip is targeted by the largest space and 
ground-based facilities for multi wavelength observations: (i) the 
DEEP and DEEP2 surveys provide optical spectroscopy to the limit 
Rab ~ 24 mag, (ii) multiwaveband optical photometry to fainter 
magnitudes is underway as part of the Canada France Hawaii Tele- 
scope Legacy Survey, (iii) deep imaging and spectroscopy, inde- 
pendent from the programs above, has been performed by Steidel et 
al. (2003) in search for Lyman break galaxies, (iv) comprehensive 
imaging with HST/ACS has recently been completed, (v) Spitzer 
mid-IR data are available, (vi) radio observations to sub-mjy levels 
have been obtained by Fomalont et al. (1991) with new much wider 
VLA observations recently completed, (vii) SCUBA has observed 
part of this field to the deep limits of the Canada-UK Deep Sub- 
millimetre survey (Webb et al. 2003). A combination of the X-ray 
observations with the mutliwavelength datasets above promises a 
breakthrough in the study of the evolution and large scale structure 
of AGNs as well as the connection between AGN activity and host 
galaxy formation. 



This paper presents the optical and X-ray spectral properties of 
the sources detected in the first 200 ks Chandra pointing observed 
as part of the Extended Groth Strip X-ray survey. This observation 
encompasses the original Groth- Westphal Strip region (Groth et al. 
1994). In addition to studying the properties of the X-ray sources 
in the context of XRB models, our purpose is to demonstrate the 
power of the full 0.5 deg 2 Extended Groth Strip Chandra survey, 
when completed, for XRB studies. Throughout this paper we adopt 
Ho = 70kms~ 1 Mpc -1 , Cl M = 0.3 and Q A =0.7. 



2 DATA 

2.1 X-ray observations 

The X-ray data used in this paper are from the A03 Chandra ob- 
servations of the original Groth Westphal Strip (GWS), which is 
part of the Extended Groth Strip (EGS) region. The total exposure 
time is about 190ks split into 3 separate integrations obtained at 
different epochs. All 3 observations were obtained with the ACIS- 
I instrument (17' x 17') with a similar roll-angle at the aimpoint, 
a = 14:17:43.6, S =52:28:41.2. A detailed description of the data 
reduction, source detection and flux estimation has been presented 
by Nandra et al. (2005). 

Briefly, standard reduction methods were applied using the 
CIAO version 3.0.1 data analysis software. After merging the in- 
dividual observations into a single event file, we constructed im- 
ages in 4 energy bands 0.5-7.0 keV (full), 0.5-2.0 keV (soft), 2.0- 
7.0 keV (hard) and 4.0-7.0 keV (ultra-hard). Source detection was 
performed using a simple but efficient method which is based on 
pre-selection of candidate sources using the WAVDETECT task of 
CIAO followed by aperture count extraction using the 90 per cent 
Point Spread Function (PSF) radius and a local background deter- 
mination to estimate the source significance. The final catalogue 
used in this paper comprises a total of 158 sources over a total sur- 
veyed area of 0.082 deg 2 to a Poisson detection probability thresh- 
old < 4 x 10~ 6 . Of these sources a total of 155, 121, 97, and 
44 are detected in the full, soft, hard and ultra-hard bands respec- 
tively. Fluxes are estimated by integrating the net counts within an 
aperture corresponding to the 70 per cent encircled energy radius 
at the position of the source. The counts in the full, soft, hard and 
ultra-hard bands are converted into fluxes in standard bands, 0.5- 
10, 0.5-2, 2-10 and 5-10keV respectively. The limiting flux in each 
of these bands is estimated 3.5 x 10" ls , 1.1 x 10~ 16 , 8.2 x 10~ 16 
and 1.4 x 10~ 15 ergs -1 cm~ 2 respectively. 

2.2 Optical photometry 

The main photometric catalogue used in this paper for the optical 
identification of the X-ray sources is the DEEP2 survey of the EGS 
that also overlaps with the original GWS field. 

The DEEP2 survey photometric data were obtained at the 
Canada-France-Hawaii Telescope (CFHT) using the 12k x 8k 
pixel CCD mosaic camera providing a 0.70 x 0.47 deg 2 field of 
view per pointing. The observations were performed in the B, 
R and I filters. The data reduction, source detection, photomet- 
ric and astrometric calibration as well as the star-galaxy separa- 
tion are described in Coil et al. (2004). The pointing that over- 
laps with the Chandra X-ray data used here is nearly complete to 
Rab ~ 24.50 mag (Bab ~ 24.75, Iab ~ 23.5 mag). This is 
shallower than the nominal limit of the full DEEP2 EGS survey 
{Rab » 24.75 mag) because of poorer seeing conditions (about 



The Deep Chandra Survey of the Groth Strip - II. 



3 



0.95 arcsec) at the time of the observations (Coil et al. 2004). The 
astrometric accuracy of the photometric catalogue is estimated to 
be 0.5 arcsec and is limited by systematic errors of the USNO-A 
catalogue used to determine the astrometric solution. 

The GWS also overlaps with the ongoing deep synoptic 
Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). This 
project uses the wide field imager MegaPrime equipped with the 
MegaCam CCD array providing a 1 x 1 deg 2 field of view. In this 
paper we use the first data obtained as part of the deep synoptic 
survey in the ugriz filters. The exposure times in each waveband 
range from 1-13 hours (depending on the filter) corresponding to 
about 2-9 per cent of the target integration at the completion of 
the project. The data reduction, source detection, photometric and 
astrometric calibration will be presented in a future paper. In brief 
the Elixir package was used for the reduction as well as the initial 
photometric and astrometric calibration, which were then refined 
using our own routines. The final astrometric uncertainty is esti- 
mated to be about 0.3 arcsec. The photometric accuracy is found 
to be better than 0.05 mag in all filters while the completeness 
limit in the AB system is r ~ 26 mag (u ~ 25.5, g ~ 26.0, 
i ~ 25.5, z ~ 25.0 mag). Although the CFHTLS deep synop- 
tic dataset reaches fainter limits that the DEEP2, we prefer to use 
the latter at present as the basic photometric catalogue because of 
its homogeneity and the well documented observational properties 
of this survey (e.g. Coil et al. 2004; Faber et al. 2005; Willmer et 
al. 2005). We nevertheless use the multi-waveband photometry of 
the CFHTLS primarily to estimate photometric redshifts but also 
to search for X-ray source optical counterparts that are fainter than 
the DEEP2 limit. 

Finally, the GWS has been targeted for deep optical imaging 
as part of a larger program searching for Lyman Break galaxies 
(Steidel et al. 2003). The observations were performed at the Kitt 
Peak 4m Mayall telescope using the Prime Focus CCD camera 
(14.2 x 14.2 deg 2 field of view) in the Un, G and K filters. Be- 
cause of the smaller field-of-view of these observations the outer 
edges of X-ray pointing do not overlap with the optical image. A 
detailed description is presented by Steidel et al. (2003). The as- 
trometry is accurate to about 0.4 arcsec and the photometric internal 
scatter is estimated to be better than 0.03 mag in all filters. These 
observations reach a limiting magnitude 5Ras ~ 26 mag, similar to 
the CFHTLS. They are used here primarily to estimate photometric 
redshifts using the Lyman break selection criteria. As discussed by 
Steidel et al. (2003), these methods are very efficient in identifying 
galaxies in narrow redshift slices in the range 1.5 < z < 3. 

We note that the ii-band filters used in the above three datasets 
are similar and therefore, there is good agreement in the estimated 
i?-band magnitudes of the same object among the different surveys. 



2.3 Optical spectroscopy 

The main source of optical spectroscopy in this study is the DEEP2 
redshift survey. This is an ongoing project that uses the DEIMOS 
spectrograph on the 10 m Keck II telescope aiming to obtain red- 
shifts for about 40000 galaxies in the range 0.7 < z < 1.5 to 
a limiting magnitude Rab = 24.1 mag. The spectra are obtained 
with a high resolution grating (12001/mm, R pb 5000) and span 
the wavelength range 6500 — 9100 A. This spectral window allows 
the identification of the Oil emission line in the redshift interval 
0-7 < z < 1.4. Outside this range the ability to measure redshifts 
and hence, the completeness of the DEEP2, drops significantly. The 
data reduction was performed using an IDL based pipeline devel- 



oped at UC-Berkeley (Cooper et al. 2006) and adapted from reduc- 
tion programs created for the SDSS. 

The GWS has also been targeted by a number of spectroscopic 
programs (Lilly 1995; Brinchmann et al. 1998; Hopkins et al. 2000; 
Voght et al. 2005) that have been compiled into a single database by 
Weiner et al. (2005 1 ). The entire EGS overlaps with the Sloan Dig- 
ital Sky Survey (SDSS) and therefore spectra for relatively bright 
galaxies and QSOs are also available (York et al. 2000). 

In addition to the above surveys Steidel et al. (2003, 2004) 
performed follow-up multi-slit spectroscopy of the GWS Lyman 
Break galaxies as well as some X-ray sources using the LRIS-B on 
the Keck telescopes. The observations used a 300 line/mm grating 
blazed at 5000A leading to a dispersion of 2.47 A/pixel, a wave- 
length range that included at least the 4000 — 7000 A regime and 
a nominal spectral resolution of about 12. 5A. The total integration 
time varied between 1.5-3 h (depending on the observing condi- 
tions), split into 1800 s sub-exposures followed by a dither of the 
telescope in the slit direction. The data were reduced using a cus- 
tom package based on IRAF scripts. A total of 10 sources in our 
sample have redshift measurements from these observations. 



3 OPTICAL IDENTIFICATION 

The main catalogue used to optically identify the GWS X-ray 
sources is the DEEP2 survey. We first search for systematic off- 
sets between the astrometric solutions of the X-ray and optical cat- 
alogues. A matching radius of 2 arcsec is adopted to include only 
secure optical identifications. We also consider X-ray sources with 
off-axis angles < 6 arcmin (total of 86) where the Chandra PSF is 
superior with a 90 per cent encircled energy radius of <4 arcsec. 
A total of 50 X-ray sources have optical identifications brighter 
than Rab < 24.5 mag. We estimate small systematic offsets of 
8RA = -0.23 and 5DEC = 0.37 arcsec between the X-ray and 
optical source positions. These were then used to align the X-ray 
source catalogue to the DEEP2 astrometric solution. 

Next we explore the positional accuracy of the X-ray centroid 
as a function of off-axis angle, 9. We match the X-ray and opti- 
cal catalogues using an ample 5 arcsec search radius to account for 
the degradation of the PSF at large off-axis angles. FigureQplots 
the positional offset in RA and DEC between the X-ray and optical 
source positions against 9. Reassuringly, the mean X-ray-optical 
offset is close to zero at all off-axis angles but the lcr rms increases 
from about 0.5 arcsec at 9 < 6 arcmin to ~ 1 arcsec at larger off- 
axis angles. We account for the degradation of the X-ray positional 
accuracy by varying the matching radius as a function of off-axis 
angle. For 9 ^ 6 arcmin we use a radius of 1.5 arcsec, correspond- 
ing to the 3u rms scatter around the mean. For 9 > 6 arcmin the 
matching radius increases to 3 arcsec, the 3<r rms positional uncer- 
tainty at these off-axis angles. 

The surface density of optical sources to the limit Rab = 
24.5 mag is large enough that a substantial fraction of chance as- 
sociations is expected within the above radii. We account for this 
effect by estimating the Poissonian probability, P, that a given op- 
tical counterpart is spurious alignment following the method of 
Downes et al. (1986). Given the surface density of objects brighter 
than m, E(< m), the expected number of candidates within r is 
/i = 7r r 2 E(< m). Assuming that source positions are Poissonian, 
the probability of at least one object brighter than m within radius 

1 http://saci.ucolick.org/verdi/public/index.html 
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r is P = 1 — exp(— /x). In practice one has to apply a cutoff in P 
to limit the optical identifications to those candidates that are least 
likely to be spurious alignments. 

The probability P however, is estimated under the assump- 
tion that the source positions are uniformly distributed within the 
surveyed area. For the real clustered distribution of optical sources 
we assess the fraction of spurious optical identifications for dif- 
ferent probability cutoffs using mock X-ray catalogues constructed 
by randomising the positions of the X-ray sources within the area 
covered by the Chandra observations. The optical identification 
method is performed on the mock catalogues using the same cri- 
teria (e.g. matching radius) as for the real sources. This procedure 
is repeated 500 times. When constructing random X-ray catalogues 
we maintain the spatial distribution of the sources due to both vi- 
gnetting and real clustering. This is accomplished by applying off- 
sets in the range 30-60 arcsec to the X-ray source positions around 
their original centroid. 

Figure|2|plots the cumulative distribution of optical identifica- 
tions for the full-band X-ray selected sample (using the matching 
radius scheme described above) as a function of probability cutoff. 
Also shown in Figure|2|is the expected number of spurious coun- 
terparts estimated as described above. This figure shows that the 
number of optically identified sources reaches a plateau at P ~ 1.5 
per cent while the spurious identification rate further increases with 
P. Based on Figure|2|we adopt a cutoff probability P < 2 per cent 
for optical identification in the case of off-axis angles < 6 arcmin. 
Because of the degradation of the positional accuracy at larger off- 
axis angles we relax the probability cutoff to P < 4 for sources 
with 8 > 6 arcmin. This is to minimise the fraction of missed opti- 
cal identifications because of the poor X-ray positions. For an op- 
tical source with Rab = 24.5 mag the probabilities P < 2 and 
P < 4 per cent correspond to maximum separations between the 
optical and X-ray centroids of about 1 and 1.5 arcsec respectively. 
Repeating the simulations above using the off-axis dependent iden- 
tification scheme we estimate a spurious fraction of about 4.5 per 
cent. The choice of P is a trade-off between maximum number of 
optical counterparts and minimum contamination rate. Similar re- 
sults and false identification rates are obtained for the soft, hard and 
ultra-hard samples. 

For X-ray sources with no optical identification to the DEEP2 
magnitude limit or outside the DEEP2 field-of-view we use the 
CFHTLS to search for fainter optical counterparts applying the 
same selection criteria described above. Moreover, a number of op- 
tically faint X-ray sources lie in the gaps between the CCDs of the 
MegaCam mosaic. For these sources we use the Steidel et al. (2003) 
deep optical imaging to search for optical identifications. Consid- 
ering sources fainter than the DEEP2 magnitude limit increases the 
spurious fraction rate by about 2 per cent. A total of 29 sources 
are identified with galaxies from the CFHTLS or the Steidel et al. 
(2003) survey. Table □ summarises the identification statistics for 
different X-ray selected subsamples. Table |2| presents the optical 
properties of the GWS X-ray sources as well as the source of opti- 
cal photometry. 




Figure 1. X-ray/optical positional offset in RA (upper panel), DEC (middle 
panel) and total angular distance (lower panel) against X-ray off-axis angle. 
The filled points represent the mean offset within different off-axis angle 
bins. The horizontal errorbar corresponds to the width of each bin while, 
the vertical errorbar is the la rms. The width of the bins varies so that each 
of them includes about 25 X-ray/optical pairs. 
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Figure 2. Number of optical identifications as a function of probability cut- 
off P for the full-band sample. The continuous histogram is for the real 
X-ray catalogue. The dashed line corresponds to the mean of 500 mock 
X-ray catalogues as described in the text. 



4 REDSHIFT ESTIMATION 

Spectroscopic redshifts are available for a total of 51 sources. These 
are classified into 3 groups on the basis of their optical spectro- 
scopic properties (primarily from DEEP1 and DEEP2) following 
methods described in Sarajedini et al. (2006): broad emission-line 
galaxies, narrow emission-line sources and systems with absorp- 



tion lines. For broad emission line AGN we adopt the criterion 
FWHM > 1200 kms -1 . We note that the low S/N ratio of some 
of the spectra and the small spectral window of the DEEP2 obser- 
vations (6000 — 9500 A) introduce some uncertainty in the clas- 
sification scheme above. In addition to the above three groups, a 
number of sources in the sample cannot be classified because the 
optical spectra were not available for visual inspection. 

For X-ray sources without spectroscopic identification we es- 
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timate photometric redshifts exploiting the multiwaveband pho- 
tometry (ugriz) of the CFHTLS. Determining photometric red- 
shifts for X-ray sources is challenging because of the significant, if 
not dominant, AGN component contributing to the optical broad- 
band colours (e.g. Babbedge et al. 2005; Kitsionas et al. 2005). Re- 
cent studies however, suggest that many moderate luminosity AGN 
(< 10 44 ergs _1 ) as well as obscured X-ray sources have optical 
continuum emission that is dominated by stellar light thus, allow- 
ing galaxy templates to be used for photometric redshifts (Barger et 
al. 2003; Gandhi et al. 2004; Georgakakis et al. 2004; Georgakakis 
et al. 2006). 

We explore this possibility using the photometric redshift code 
of Gwyn (2001) based on a standard \ 2 minimisation method. The 
galaxy templates are based on those of Coleman, Wu & Widman 
(1980), providing SEDs for 4 main spectral galaxy types (E/S0, 
Sbc, Scd, Im), extended in the UV and IR wavelength regions us- 
ing the GISSEL98 code (e.g. Bruzual & Chariot 1993). These are 
supplemented with the starburst SB2 and SB3 spectra from Kinney 
et al. (1996). Having only a small number of SEDs can cause alias- 
ing in photometric redshifts. Therefore, a new template set has been 
created by smoothly interpolating between each of the six original 
spectra. This results in a total 51 SEDs from ellipticals (spectral 
classification 0) to extreme starbursts (spectral classification 1). 

Figure[3]compares the photometric and spectroscopic redshift 
estimates for the sources with available spectroscopic observations. 
With the exception of broad-line QSOs, for which galaxy templates 
are inappropriate for photometric redshift estimates, there is fair 
agreement between z 3pec and z p hot with an rms scatter, after ex- 
cluding broad-line AGN, i £ ( ^ ~*^ cc ) = 0.08, where 
N is the total number of sources. This is further demonstrated in 
Figure |4| plotting Sz — (z p hot — z sp ec) against B — I colour. 
Sources bluer than B — I ~ 1.5 are dominated by broad-line 
QSOs and show significant scatter in their photometric redshift de- 
termination (e.g. Barger et al. 2003; Kitsionas et al. 2005). For sys- 
tems redder than this limit however, the photometric redshifts are 
more reliable. To avoid erroneous redshift estimates in the analy- 
sis that follows we use photometric redshifts only for sources with 
B - I > 1.5 mag. 

In addition to the standard photometric redshift estimation 
above, we exploit the Lyman Break galaxy selection available for 
the GWS (Steidel et al. 2003, 2004) to determine the redshift of a 
small number of X-ray sources. This method has been shown to be 
very efficient in identifying galaxies in well defined narrow redshift 
slices within the range 1.5 < z < 3 with a low interloper rate. Here 
we use the BM, BX, C, D, M and MD Lyman Break galaxy selec- 
tion criteria fully described in Steidel et al. (2003, 2004), which 
correspond respectively to redshifts 1.70 ± 0.34, 2.20 ± 0.34, 
3.09 ± 0.22, 2.93 ± 0.26, 3.15 ± 0.24 and 2.79 ± 0.27. In our 
sample there are 7 X-ray sources that fullfill one of the above se- 
lection criteria: 2 BX, 2 MD, 1 C, 1 D and 1 M. 

The optical spectroscopic and photometric redshift informa- 
tion for the GWS sample is presented in Table [2] For the full and 
hard band samples Figure|5]presents the optical magnitude distri- 
bution of sources with spectroscopic, photometric or no redshift 
information as well as sources without optical identification. The 
total number of sources in these groups for different X-ray selected 
samples is also shown in Table Q The no-redshift class involves 
sources that are either too faint to estimate photometric redshifts, 
do not have CFHTLS data (e.g. CCD gaps) or have counterparts 
with B — I < 1.5 mag and therefore unreliable redshift determina- 
tion. 
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Figure 3. Photometric against spectroscopic redshift estimates for the 
X-ray sources with available spectroscopic observations. Circles are for 
sources with absorption-line spectra, squares correspond to systems with 
narrow emission-line spectra and crosses are broad line AGNs. Triangles 
are sources with no classification. 
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Figure 4. 8z = (z p hot ~ Zs P ec) against DEEP2 B — I colour. The symbols 
are the same as in Figurel3l 

5 X-RAY SPECTRA 

For the X-ray spectral analysis we use the XSPEC vl 1.3.1 package. 
The X-ray counts of each source are extracted using the 95 per cent 
encircled energy radius (1.5 keV) at the position of the source. The 
background is estimated using an annulus centered on the source 
with inner aperture size 1.5 times larger than the 95 per cent en- 
circled energy radius and outer aperture 100 pixel greater. We fit 
the data adopting a power-law model absorbed by both an intrinsic 
column density at the redshift of the source and a Galactic column 
at z = fixed to Nh = L3 x 10 20 cm -2 , appropriate for the EGS 
(WABS*ZWABS*POW). For the absorption we adopt the Wiscon- 
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X-ray total optical spectro-z photo-z optical ID, no-ID 

sub-sample number IDs only only no-z 



Total 


155 


128 


51 


36 


42 


27 


Soft 


121 


102 


42 


27 


33 


19 


Hard 


97 


83 


37 


24 


22 


14 


Ultra-hard 


44 


40 


23 


10 


7 


4 



Table 1. Optical and spectroscopic identification statistics. The columns are: (1): X-ray sub-sample; (2) total number of X-ray sources; (3) number of optical 
counterparts; (4) number of spectroscopic identifications; (5) number of sources with photometric redshift determination only (i.e. not spectroscopic z); (6) 
number of sources with optical counterparts for which photometric redshift estimations was not possible; (7) number of blank fields. 
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Figure 5. Rab optical magnitude distribution for different subsamples 
of X-ray sources. The four panels from bottom to top plot respectively 
histograms for sources with spectroscopic redshifts, photometric redshifts 
only, without any redshift determination (but with an optical ID) and with 
no optical identification. Within each panel the continuous line corresponds 
to the full-band sample and the dashed line (offset to the right by 0. 1 mag 
for clarity) represents the hard sample. 



sin cross-sections (Morrison and McCammon 1983). For sources 
without spectroscopic or photometric redshift estimates we assume 
z = 1.5 to estimate the intrinsic Nh- Adopting a different mean z 
in the range 1 — 2 for these systems however, does not significantly 
modify our results and conclusions. 

In the case of sources with small number of net counts (< 200) 
we use the C-statistic technique (Cash 1979) specifically devel- 
oped to extract information from low signal-to-noise ratio spectra. 
The data are grouped to have at least one count per bin. We at- 
tempt to constrain the intrinsic Nh by fixing the power-law index 
to r = 1.8. This value of Y is selected to be inbetween the mean 
spectral index of radio loud (r = 1.6; Reeves & Turner 2000; 
Gambill 2003) and radio quiet AGNs (T » 1.9; Nandra & Pounds 
1994; Reeves & Turner 2000). Adopting a single fixed F for the 
spectral analysis is an approximation. We caution that our analy- 
sis will overestimate the Nh for sources with spectra intrinsically 
flatter than V = 1.8. 

For sources with sufficient number of counts (> 200) we per- 
form standard x 2 spectral fitting. The data were grouped to have 
a minimum of 20 counts per bin to ensure that Gaussian statistics 



apply. For the \ 2 analysis we require that the source spectrum has 
at least 10 spectral bins. The WABS*ZWABS*POW model provides 
acceptable fits (i.e. reduced \ 2 ~ 1) for all sources. The parame- 
ters estimated from the C-statistic and the % 2 analysis are consistent 
within the errors. 

For both the % 2 an d the C-statistic analysis the fit was per- 
formed in the 0.5-8 keV energy range where the sensitivity of the 
Chandra is the highest. The estimated errors correspond to the 90 
per cent confidence level. The results of the X-ray spectral analysis 
are presented in Table|2| 



6 THE MODEL 

In this section we describe the model we use to interpret the opti- 
cal and X-ray properties of the GWS X-ray sources in the context 
of AGN evolution scenarios and different parameterisations for the 
intrinsic Nh distribution. Modeling of the data requires certain as- 
sumptions about the X-ray spectra of AGN, their luminosity func- 
tion and its evolution with redshift as well as the relative fraction 
of obscured and unobscured systems. 

We model the X-ray spectra of AGN adopting for simplic- 
ity an absorbed power-law spectral energy distribution with fixed 
exponent T — 1.8 and photoelectric absorption cross sections as 
described by Morrsion & McCammon (1983) for solar metallicity. 

For the X-ray luminosity function (XLF) of AGNs and its evo- 
lution with redshift we use the two different parameterisations pre- 
sented by Miyaji et al. (2000) and Ueda et al. (2003). 

Miyaji et al. (2000) combined deep pencil-beam and shal- 
low wide-area ROSAT surveys to estimate the XLF of unobscured 
AGNs in the rest-frame 0.5-2 keV energy band. We adopt the lu- 
minosity dependent density evolution parameterisation of the XLF 
proposed by Miyaji et al. (their model LDDE1). These authors ar- 
gue that this model provides a better description of the observations 
compared to the pure density or luminosity evolution. In this picture 
AGNs evolve differentially with more luminous systems evolving 
faster than less luminous ones. Such a trend has also been proposed 
for the optical luminosity function of QSOs (e.g. Wisotzki 1998). 

Ueda et al. (2003) estimated the AGN XLF in the rest-frame 2- 
lOkeV energy range using a combination of hard-band (> 2keV) 
surveys conducted with HEAO-1, ASCA and Chandra missions. 
Here we adopt the luminosity dependent density evolution of the 
luminosity function, which according to Ueda et al. provides a bet- 
ter fit to the data. In this parameterisation the cutoff redshift, after 
which the evolution of AGN stops, increases with luminosity. We 
note that this is different from the Miyaji et al. (2000) luminosity 
dependent density evolution, where it is the rate of evolution that 
changes with luminosity but not the cutoff redshift. 

For the AGN Nh distribution, /(Nh), we experiment with 
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different parameterisations. The first model adopted here is the one 
estimated by Ueda et al. (2003) on the basis of observational data 
and for column densities in the range 10 zu < N H < cm 
The interesting feature of their functional form is that the frac- 
tion of obscured AGNs drops with increasing luminosity. In that 
sense the Ueda et al. (2003) /(Nh) does not strictly follow the 
unified model prescription (Antonucci 1993) where the only pa- 
rameter determining the obscuring column density is the viewing 
angle to the observer. We also consider a set of models where 
f(Nn) is a step function with a fixed ratio, 1Z, between ab- 
sorbed (N H > 10 22 cm -2 ) and unabsorbed (N H < 10 22 cm" 2 ) 
sources. We further assume that obscured and unobscured AGNs 
are distributed uniformly in the range 10 20 < N H < 10 22 and 
10 22 < N H < 10 26 cm respectively. For the obscured sys- 
tems, the above assumption is in fair agreement with the Nh dis- 
tribution of Seyfert-2s estimated by Risaliti, Maiolino & Salvati 
(1999). These authors find that about 75 per cent of their sample has 
N H > 10 23 cm" 2 and at least 25 per cent has N H > 10 25 cm" 2 . 
This set of models are consistent with the unified scheme with the 
ratio 1Z related to the opening angle of the torus. Similar models but 
with a smooth transition between obscured and unobscured sources 
are presented by Treister et al. (2004). In this paper we use models 
with 1Z = 1, 2, 3 and 4. We note that 1Z — 4 is the locally estimated 
value for the fraction of obscured AGN (Maiolino & Rieke 1995). 



7 RESULTS 

In this section we compare the observed redshift and column den- 
sity distributions of the GWS sample against the predictions of the 
models above. To minimise incompleteness uncertainties due to ei- 
ther optically unidentified sources or systems without redshift de- 
termination (spectroscopic or photometric) we perform the com- 
parison for the hard-band X-ray sources. The 2-10keV sample has 
indeed, sufficient number of sources to avoid poor statistics while, 
about 37 per cent of them (36 out of 97) are either blank fields 
or do not have a redshift estimate. Selection at the 2-10keV band 
also provides samples that are less sensitive to obscuration and is 
therefore best suited for studies on the intrinsic fraction of obscured 
AGN. Finally, many groups have published results for the 2-10 keV 
spectral band and therefore, choosing this energy range for the anal- 
ysis facilitates comparison of our survey with previous samples. For 
the comparison with the models above, the AGN XLF is integrated 
in the redshift interval z = — 5 for unobscured luminosities in 
the range L x (2 - lOkeV) = 10 42 - 10 46 ergs" 1 . For each lu- 
minosity and redshift interval the distrubition of AGNs to different 
columns is described by the /(Nh) models discussed in the previ- 
ous section. The predicted number density of objects in each Lx, 
Nh and z bins are then folded through the sensitivity map of the 
EGS to estimate the number of AGNs to the flux limit of the survey. 



7.1 The redshift distribution 

Figure |6| shows the redshift distribution of the hard-band sample. 
For sources without spectroscopic redshifts we use the photometric 
redshift probability density distribution, instead of the primary so- 
lution only, to construct the histogram in Figure [6] This approach 
guarantees that some of the uncertainties involved in the determi- 
nation of photometric redshifts are factored into our analysis. For 
sources assigned photometric redshifts based on the Lyman break 
galaxy selection (Steidel et al. 2003, 2004) we assume a Gaussian 



probability density distribution with a mean and a standard devi- 
ation appropriate for the selection criteria that each source fulfills 
(see section|4j- 

A total of 36 hard X-ray selected sources (37 per cent of the 
sample) do not have spectroscopic or photometric redshift determi- 
nation. These are shown with the hatched histogram in Figure [6] 
Their optical magnitude distribution is presented in the two upper 
panels of Figure [5] Fourteen of these 36 sources are blank fields 
and are most likely associated with z > 1 systems. The remaining 
22 sources have a distribution that is skewed to fainter magnitudes 
in Figure [5] compared to spectroscopically identified systems but 
similar to that of X-ray sources with photometric redshift determi- 
nation. However, 17 of these 22 sources have B — I < 1.5 mag and 
are likely to be associated with high-z QSOs. Nevertheless, unless 
all spectroscopically unidentified sources are clustered in a narrow 
redshift slice, we do not expect them to drastically modify the po- 
sition of the peak of the distribution in Figure[6] 

Also shown in Figure [6] are the predictions of the two model 
XLFs presented in the previous section for the Ueda et al. (2003) 
/(Nh)- The adopted Nh distribution has only minor effects on 
the resulting redshift distribution and does not affect any of our 
conclusions. In Figure|6|the Miyaji et al. (2000) prediction peaks 
at z ~ 1.5, higher than the observations. On the contrary, the Ueda 
et al. (2003) XLF produces a redshift distribution with a peak and 
overall shape in broad agreement with the data. There is however, 
a larger fraction of z w 1 sources compared to the model pre- 
diction, suggesting a comsic variance spike. This is more clearly 
demonstrated in the inlet plot of Figure [6] which uses a narrower 
logarithmic redshift bin of 0.05. The full EGS sample will have a 
sufficiently wide FOV (0.5 deg 2 ) to address this issue. The Ueda 
et al. luminosity function also predicts a larger number of high- 
z systems compared to the observations. It is possible that some 
the X-ray sources without redshift determination will populate this 
high-z tail. We attempt to quantify the agreement between the ob- 
served and model distributions in Figure|6| in the optimal case that 
the spectroscopically unidentified sources are distributed to redshift 
bins in such a way that the difference between the observed and 
model N(z) is minimal. In the case of the Ueda et al. (2003) XLF 
we estimate a x 2 -test probability that the two distributions (model 
and observations) are drawn from the same parent population of 
about 99 per cent. For the Miyaji et al. (2000) model this exercise 
gives a probability of < 1 per cent. 

7.2 The column density distribution 

Figure Q presents the Nh distribution of the hard X-ray selected 
sample. The fraction of obscured AGN (Nh > 10 22 cm" 2 ) in this 
figure is estimated 48 ± 9 per cent. For sources without spectro- 
scopic or photometric redshift estimates we assume z — 1.5 to 
estimate the intrinsic Nh ■ Adopting a different mean z in the range 
1 — 2 for these systems however, does not significantly modify the 
derived Nh distribution. Also, the photometric redshift uncertain- 
ties have little impact on the observed distribution in Figure As- 
suming a redshift dependence of the rest-frame column density of 
the form N H oc (1 + z) 2 65 (e.g. Barger et al. 2003) and a photo- 
z rms scatter 5z/(l + z) = 0.08 (see section 4), we estimate 
51ogNH = 0.09, which is much smaller than the small photon 
statistics uncertainty. 

Also, this plot is constructed using the Nh probability density 
distribution for each source instead of the best-fit solution. The ad- 
vantage of this approach is that the column density uncertainties, 
due to small number of photons in the X-ray spectra, are taken 
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Figure 6. Redshift distribution of the hard X-ray selected sample. The 
hatched region shows the number of sources without spectroscopic iden- 
tification. The dotted and dashed lines are the predictions of the Ueda et al. 
(2003) and Miyaji et al. (2000) XLFs respectively. The inlet plot shows the 
same redshift distribution only for sources in the range < z < 2 using a 
narrower bin size to search for cosmic structures within the surveyed area. 



into account in the analysis. Indeed, the C-statistic, used in the 
few counts limit, provides confidence limits and probability den- 
sity distributions for the fitted parameters. In particular the AC is 
distributed as \ 2 with v degrees of freedom (where v is the num- 
ber of fitted parameters) and hence, the same methods used to es- 
timate confidence intervals for the model parameters in the case of 
the \ 2 analysis also apply to the C-statistic for spectra with few 
counts (Cash 1979). We confirm this by performing simulations of 
power-law spectra with fixed T = 1.8 and different levels of X- 
ray obscuration. We verify that even to the limit of less than « 10 
photons the Nh probability density distribution estimated from the 
C-statistic is in agreement with that derived from the simulations. 

FigureQpresents the predictions of the different f(Nu) mod- 
els discussed in the previous section combined with the Ueda et al. 
(2003) XLF. Comparison with the data suggests that the TZ = 2 
model is in broad agreement with the observations. The x 2 -test 
probability that the two distributions are drawn from the same par- 
ent population is about 5 1 per cent. Higher or lower values of TZ 
provide poorer fits and fail to predict the observed distribution. Us- 
ing the x 2 statistical test we estimate probabilities of 12 per cent for 
the TZ = 3 model and < 1 per cent for the TZ = 1 and 4 models. For 
the Ueda et al. (2003) luminosity dependent f(Nn) we estimate a 
X 2 -test probability of about 71 per cent, somewhat better than the 
TZ = 2 model, the best of the step function Nh distributions with 
fixed ratio TZ. We note that the adopted XLF has little effect on the 
model Nh distributions shown in FigureQand therefore, does not 
alter our main conclusions. 

We note however, that our sample may comprise a number 
of Compton thick AGN (N H > 10 24 cm" 2 ) where the direct X- 
ray emission is completely blocked from view and the spectrum in 
the Chandra energy band is a pure reflection continuum. Fitting a 
single absorbed power-law to these sources is clearly not appro- 
priate and will produce erroneous Nh estimates. We attempt to 
quantify this effect by simulating reflection dominated spectra and 



Figure 7. Nh distribution of the hard X-ray selected sample. In both panels 
the bold continuous line is the observed distribution. For clarity the compar- 
ison with the model predictions is split into two panels. Upper panel: com- 
parison of the observations with the 1Z=2, 3, and 4 models described in the 
text corresponding to the dash-dotted, dashed and dotted lines respectively. 
Lower panel: comparison of the observations with the 1Z=l model (dotted 
line) and the luminosity dependent Nh distribution derived by Ueda et al. 
(2003; dashed line). All models are for the Ueda et al. (2003) XLF. 



then fitting them with WABS*POW XSPEC models as described in 
section|5| For this exercise we use the Compton reflection models 
of Magdziarz & Zdziarski (1995) as implemented in the PEXRAV 
spectral energy distribution of XSPEC. We assume a solid angle 
of 27r, solar abundance for all elements and an average inclina- 
tion relative to the line of sight cosi = 0.45. Only the reflection 
component was used, i.e. no direct radiation. We also add a FeKa 
iron line assuming a Gaussian profile with width a — 100 eV, 
similar to the instrumental FWHM of the ACIS-I and rest frame 
equivalent width of 1 keV appropriate for heavily obscured AGNs 
(N H > 10 24 cm~ 2 ; e.g. George & Fabian 1991). For the simu- 
lations we adopt a redshift z = 1.5 and fix the normalisation so 
that the spectrum has a flux of about 5 x 10~ 15 ergs -1 cm -2 in 
the 2-10keV band. The simulated spectra are fit with an absorbed 
power-law as described in section [5] to estimate the Nh our X- 
ray spectral analysis produces for this type of sources. The simula- 
tions give a narrow distribution for the estimated Nh in the range 
10 23 - 10 24 cm . We therefore underestimate the column density 
of these systems and it is likely that some of the sources in the range 
Nh ~ 10 23 — 10 24 cm -2 in FigureQshould be moved to higher 
Nh values. We note however, that even if some of the sources in 
this column density range are reflection dominated Compton thick 
AGN, we do not expect this to modify our conclusions about the 
agreement between the data and different model Nh distributions. 
Unfortunately, the small number of counts in most of the obscured 
sources in our sample does not allow us to identify Compton thick 
AGN candidates dominated by reflection emission. 



8 DISCUSSION 

In this paper we explore the redshift and the intrinsic Nh dis- 
tributions of AGN using a deep 200 ks Chandra pointing in the 
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Groth- Westphal Strip region. This is the first of a total of 8 obser- 
vations that are currently underway as part of a deep wide angle 
(0.5 deg 2 ) X-ray survey in the Extended Groth Strip. A wealth of 
optical photometric and spectroscopic data are available in this field 
(e.g. DEEP2, CFHTLS) providing optical identifications as well as 
spectroscopic and photometric redshift estimates for the X-ray pop- 
ulation. 

The advantage of this dataset is that the detected sources are 
responsible for a sizable fraction of the XRB (about 70 per cent 
in the 2-10keV band) and have mean X-ray spectral properties 
consistent with the X-ray background (T « 1.4; Nandra et al. 

2005) . Nevertheless they are, on average, brighter than the ex- 
tremely faint X-ray population identified in the ultra-deep Chan- 
dra surveys (Alexander et al. 2003) facilitating follow-up multi- 
wavelength studies. For example, the hard X-ray selected sample 
comprises 97 sources of which 74 have R < 24.5 mag (76 per 
cent) and therefore are accessible for optical spectroscopy using 
10-m class telescopes. For comparison, in the CDF-N a total of 332 
sources are detected in the 2-8 keV spectral band of which 203 (60 
per cent) have R < 24.5 mag and 162 (50 per cent) have optical 
spectroscopy available (Barger et al. 2003). 

For comparison with the models we focus on the hard X-ray 
selected sample. This combines sufficient number of sources (97) 
for statistical reliability and high optical identification rate min- 
imising uncertainties due to optically faint X-ray sources. A total 
of 14 systems in this sample (14 per cent) are blank fields. Many of 
them also have log fx/fopt > 1. This is shown in Figure|S|which 
plots i?As-band magnitude against 2-10 keV X-ray flux. High X- 
ray-to-optical flux ratio systems are suggested to comprise a large 
fraction of high- 2 heavily obscured type-II QSOs (Mignoli et al. 
2004; Civano, Comastri & Brusa 2005). Figure[9|plots the N H dis- 
tribution of these sources assuming a mean redshift z = 1.5. The 
optically unidentified sources in the hard sample are skewed to- 
ward high column densities compared with identified sources, with 
a median of about 5 x 10 22 cm -2 . 

The redshift distribution of the GWS hard-band sample, us- 
ing both spectroscopic and photometric redshifts, shows a peak at 
2 ~ 1. Although about 2/5 of the X-ray sources do not have redshift 
determination, we argue that these systems are unlikely to modify 
the shape of the N(z). The observed redshift distribution in Fig- 
ure |6|is in agreement with previous deep surveys in the hard-band 
that also find that the X-ray population peaks at z ~ 1 (Fiore et al. 
2003; Georgantopoulos et al. 2004; Barger et al. 2005; Treister et 
al. 2005). This is at odds with the recent versions of the population 
synthesis models that successfully reproduce the spectral proper- 
ties of the XRB but predict a peak at higher redshifts z ~ 1.5 
(Comastri et al. 1995; Gilli et al. 2001). As demonstrated in Fig- 
ure IH the origin of this discrepancy is that the models above adopt 
the XLF derived for soft-band selected powerful QSO that peak at 
2 ~ 1.5 — 2 (e.g Miyaji et al. 2000). More recent studies however, 
that combine X-ray selected AGN samples over a wider luminosity 
range, suggest a more complex evolutionary history that strongly 
depends on Lx '■ more powerful systems evolve up to z ~ 2 while 
less luminous sources peak at lower redshift, z ~ 1 (Ueda et al. 
2003; Barger et al. 2005; Hasinger, Miyaji & Schmidt 2005). As 
shown in Figure [6] such a luminosity dependent density evolution 
can successfully reproduce the redshift distribution of our sample. 

The above complex evolutionary pattern warrants some phys- 
ical interpretation. It is possible that less luminous systems, sug- 
gested to comprise a higher fraction of type-II AGN (Ueda et al. 
2003; Barger et al. 2005; Akylas, Georgantopoulos & Georgakakis 

2006) , are linked to starburst activity with a peak at z « 1, while 
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Figure 8. R^s-band magnitude against 2-10keV flux for the hard-band 
selected sample. The lines log fx / fopt = ±1 delineate the region of the 
parameter space occupied by powerful unobscured AGNs and are estimated 
from the relation log f x /fopt = log fx (2 - 10 keV) + 0.4 R AB + 5.46. 
A cross on top of a symbol is for sources with broad-line optical spectra. 
Open squares and open circles on top of a dot correspond to sources with 
narrow emission-line and absorption optical spectra respectively. Triangles 
represent upper limits in optical magnitude for sources without optical iden- 
tifications. 
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Figure 9. Nh distribution of the hard X-ray selected sources with no opti- 
cal identification. 



more powerful sources follow the QSO evolution peaking at z ~ 2. 
In this direction, the large fraction of hard (5-10keV) X-ray se- 
lected sources with mid-infrared counterparts (Fadda et al. 2002) 
has motivated models where obscured X-ray sources are tied to the 
infrared-luminous population with an evolution that steeply rises 
to 2 « 1 and levels off at higher-2 (Franceschini, Braito & Fadda 
2002; Gandhi & Fabian 2003). These models produce a redshift 
distribution for the faint X-ray population with a peak at z m 1 and 
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Figure 10. fraction of sources with Nh > 10 22 cm -2 as a function of 2- 
lOkeV flux. Open triangles: Chandra Deep Field South adapted from Aky- 
las et al. (2006); Open circles: the XMM-Newton survey described by Aky- 
las et al. (2006); Filled circles: this study. The thin-line curves from bottom 
to top correspond to models with TZ = 1 — 4. The bold curve is the model 
that uses the Ueda et al. (2003) Nh distribution. 



uncertainties. Also plotted in this figure are the predictions of the 
different /(Nh) models described in section[6|adopting the Ueda 
et al. (2003) XLF. We note that the adopted XLF has little effect on 
the model curves plotted in Figure fTol At bright fluxes the observa- 
tions are in better agreement with TZ ~ 1, while at fainter fluxes the 
data are progressively more consistent with higher TZ models. This 
may suggest a luminosity dependent Nh distribution similar to that 
proposed by Ueda et al. (2003). This model is indeed in fair agree- 
ment with our data at intermediate fluxes but is not as successful 
at both the faint and the bright end of Figure fTOl although some- 
what better than a simple fixed- TZ Nh distribution. This may sug- 
gest a steeper luminosity dependence of the obscured AGN fraction 
compared to the Ueda et al. (2003) parametrisation: i.e. » 0.8 at 
Lx ~ 10 42 ergs -1 decreasing to about 0.2 at Lx ~ 10 45 ergs -1 
compared to ~ 0.6 and 0.3 respectively for the Ueda et al. (2003) 
model. We note that a number of recent studies, using ultra-deep 
and/or shallow wide-angle samples, also argue in favor of a lumi- 
nosity and/or redshift dependent Nh distribution to explain the ob- 
served properties of the X-ray population (Treister & Urry 2005; 
La Franca et al. 2005; Akylas et al. 2006). 

A luminosity dependent fraction of obscured AGN can be un- 
derstood in terms of modified unification schemes. In these models 
the inner radius and/or the geometric height of the torus vary with 
the power of the central engine because of dust evaporation and 
radiation pressure (e.g. Lawrence 1991; Simpson 2005). 



also provide acceptable fits to both the spectral shape of the XRB 
and the source counts. We note that a possible association of X- 
ray obscuration and star-formation activity has also been proposed 
to interpret the sub-mjy radio properties of X-ray selected AGNs 
(Bauer et al. 2002; Georgakakis et al. 2004). The models above 
however, require revision of the basic assumption of the unification 
model that requires that type-I and II sources evolve in lockstep, 
since they are drawn from the same parent population. 

A key issue in the study of the origin of the XRB is the in- 
trinsic Nh distribution of AGN. Even moderate amounts of gas 
(Nh ~ 10 22 cm -2 ) have a strong effect on the X-ray emission be- 
low about 2 keV and therefore this issue is better addressed by ob- 
servations at harder energies. In this respect, X-ray surveys in the 
2-10keV band at relatively bright fluxes (> 10 -14 ergs -1 cm -2 ) 
suggest that the observed fraction of AGN with Nh > 10 22 cm -2 
is inconsistent with the locally determined fraction of 4: 1 (Maiolino 
& Reike 1995) and closer to a ratio of about 1:1 (Fiore et al. 
2003; Perola et al. 2004; Georgantopoulos et al. 2004; Akylas et 
al. 2006). Therefore, XRB models that adopt a fixed 4: 1 fraction 
of obscured AGNs independent of redshift or luminosity have a 
problem reproducing the observations at bright fluxes (e.g. Treis- 
ter & Urry 2005; La Franca et al. 2005; Akylas et al. 2006). Sur- 
prisingly, deeper X-ray surveys show a different picture, suggest- 
ing an abrupt increase of the fraction of obscured AGNs at faint 
fluxes (< 10 -15 ergs -1 cm -2 ; Perola et al. 2004; Treister et al. 
2005; Akylas et al. 2006). This is demonstrated in Figure [Tol plot- 
ting X-ray flux against fraction of AGN with Nh > 10 22 cm -2 
for both the bright XMM-Newton survey of Akylas et al. (2006) 
and the CDF-South. We note that for the latter field the Nh of 
each source is estimated using spectral fitting instead of hardness 
ratio. A full description of the spectral analysis in the CDF-South 
is presented by Akylas et al. (2006). In the same figure, the re- 
sults from the EGS survey, split into two independent flux bins, are 
in broad agreement with the Akylas et al. (2006) study within the 



9 CONCLUDING REMARKS 

This paper presents first results from an ongoing deep (200 ks per 
pointing) X-ray Chandra survey in the Extended Groth Strip re- 
gion. Data from the first of a total of 8 pointings in this field are 
used here. We analyse and discuss the optical and X-ray properties 
of the sample in the context of XRB population synthesis models. 

We first construct the photometric and spectroscopic redshift 
distribution of the 2-10keV selected sample to the limit ~ 8 x 
10 -15 erg s -1 cm -2 which has a peak at z ~ 1 in agreement with 
previous studies. Although there is about 40 per cent redshift in- 
completeness we find that luminosity dependent density evolution, 
where lower luminosity systems peak at lower redshifts, is in fair 
agreement with the observations. Luminosity functions that assume 
evolution out to z m 2 for the entire AGN population (both type-I 
and II), similar to that of unobscured broad-line QSOs do not fit the 
data. The luminosity dependent evolution supported by our data is 
consistent with scenarios suggesting that lower-luminosity systems 
(dominated by obscured AGNs) are associated with star-formation 
activity and peak at z ~ 1 while, more powerful QSOs evolve out 
to higher-z. 

We also explore the Nh distribution of the sample which is 
consistent with either a fixed obscured AGN fraction of TZ ~ 2 
or the luminosity dependent Nh distribution proposed by Ueda et 
al. (2003), where less luminous systems comprise a higher fraction 
of type-II AGNs. We also argue that such luminosity dependent 
parametrisation of the Nh distribution is essential to account for 
the fraction of obscured AGN observed in different samples over a 
wide range of flux limits. 

The X-ray survey of the Extended Groth Strip region, when 
completed, will cover a total of about 0.5 cleg 2 at flux limits simi- 
lar to those presented here. The final sample will be about 8 times 
larger, comprising a total of over 1000 sources. The wealth of 
follow-up optical photometric and spectroscopic data available in 
this field (DEEP2, CFHTLS) will provide optical identifications as 
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well as spectroscopic and photometric redshift estimates for a large 
fraction of the X-ray population, for example close to 90 per cent 
for the hard-band. The use of Spitzer mid-infrared data will further 
increase the number of X-ray identifications and also has the poten- 
tial to refine the photometric redshift estimates, particularly for the 
optically faint subsample. Such a high identification rate combined 
with photometric/spectroscopic redshifts and the complementary 
multi-wavelength data (infrared, radio, sub-mm) available for the 
EGS promise a major step forward in our understanding of the na- 
ture and the evolution of the AGN populations that make up the 
XRB. 
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Table 2. Chandra GWS X-ray/optical catalogue. Col.(l): Source catalogue number; Col.(2): X-ray right ascension in J2000; Col.(3): X-ray declination in 
J2000; Col.(4): X-ray/optical centroid offset in arcsec; Col.(5): percent probability the optical counterpart is spurious alignment; Col.(6): Rab magnitude 
and source of optical identification 1=DEEP2, 2=Steidel et al. (2003), 3=CFHTLS; Col.(7): spectroscopic redshift and source catalogue 1=DEEP2, 2=DEEP 
(Weiner et al 2005), 3: CFRS (Lilly et al. 1995), 4: Steidel et al. (2003), 5: SDSS; Col.(8): optical spectroscopic classification, NL=narrow emission lines, 
BL=broad emission lines, AB=absorption lines, UNCL=no classification available, STAR=Galactic star; Col. (9): photometric redshift; Col. (10): Nh in units 
of 10 22 cm -2 ; Col. (11): unobscured 0.5-10 keV X-ray luminosity in ergs - 1 . Col. (12); Flags fshu=source detected at < 4 X 10~ 6 probability in this band, 
where the bands are f=full, s=soft, h=hard, u=ultrahard. 
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Table 2 - continued 
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Table 2 - continued 



Cat 


a x 


Sx 


Sox 


P 


Rab 


%spec 


class 


Zphot 




L x (0.5 - lOkcV) 


Flags 


No. 


(J2000) 


(J2000) 


(arcsec) 


(per cent) 


(mag) 








(10^ 2 cm 2 ) 


(10 43 ergs -1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


clOl 


14 17 56.87 


+52 31 24.49 


0.09 


0.01 


23.98 1 






740+ ' 09 


oi+°- 43 

u ' u -0.01 


1.01 


fshu 


cl02 


14 17 56.92 


+52 31 18.47 


0.00 


0.44 


25. 10 3 








2.23 +43 ^I 5 


0.53 


f 


cl03 


14 17 57.12 


+52 26 30.98 


0.68 


< 0.01 


24.45 2 








7.74+ 3 -°° 


14.69 


fshu 


cl04 


14 17 57.48 


+52 31 06.92 


0.00 


0.19 


24. 89 3 


3.026 4 


UNCL 


3.150+°? 4 

-0.24 


20.02 +2 °- t ;5 
15.41 


14.62 


fshu 


cl05 


14 17 57.51 


+52 25 46.45 


0.15 


0.02 


23. 34 1 


0.995 1 


NL 




J — 1.29 


1.51 


fsh 


cl06 


14 17 58.17 


+52 31 33.63 


0.00 


1.64 


23. 56 3 








0.12+ 4 " 
v — 0.12 


0.42 


f 


cl07 


14 17 58.20 


+52 21 53.24 


0.00 


0.32 


22.89 3 






0.940+°^ 

— U.Uo 


0.65+179 

— U.DD 


1.80 


fsh 


cl08 


14 17 58.97 


+52 31 38.89 


0.20 


0.02 


22.71 1 


0.644 4 


BL 




0.29+?- 3 * 
v '^" / — 0.29 


2.17 


fshu 


cl09 


14 17 59.32 


+52 24 20.30 


0.22 


0.07 


23.78 1 






0.320+ ,'??, 

— U. Iz 


< 0.41 


0.04 


fs 


cllO 


14 18 00.08 


+52 22 23.26 


0.24 


0.20 


24.63 1 






1.170+ ,-; 2 

' —0.13 


4.29+ 4 - 74 

^—4.29 


11.05 


fshu 


clll 


14 18 00.41 


+52 28 22.22 


0.61 


0.55 


23.73 1 








3.11+r?'? 

J -3.11 


2.33 


fh 


cll2 


14 18 00.43 


+52 36 10.11 


0.00 


1.16 


26.23 3 








< 1.56 


5.55 


fsh 


cll3 


14 18 01.15 


+52 29 41.78 


0.06 


< 0.01 


23.04 1 


2.907 4 


BL 


2.930 + n 2^ 
Urjyj — 0.2b 


< 6.27 


5.42 


fs 


cll4 


14 18 01.37 


+52 31 50.77 






<26.00 






1 090 + SL 

UUKJ —0.09 


10 21+i 4 ;° 4 


1.20 


fh 


cll5 


14 18 01.67 


+52 28 00.66 


0.02 


< 0.01 


24.67 1 








u - 00 — 0.88 


0.74 


fs 


cll6 


14 18 02.00 


+52 35 14.53 


0.32 


< 0.01 


19.93 1 


1.497 1 


BL 


2.879+°' 2 ? 


< 2.57 


111.24 


fshu 


cll7 


14 18 02.41 


+52 21 32.36 


0.00 


0.01 


18.86 3 






230+° : °? 


< 0.27 


0.05 


fsh 


cll8 


14 18 02.93 


+52 35 47.11 


0.00 


0.26 


24.97 3 








59+3.1 2 

u u — 0.59 


8.09 


fshu 


cll9 


14 18 04.55 


+52 36 33.15 


0.43 


0.28 


23.57 1 






0.810+0- 29 

— U . U4 


— U.b4 


4.81 


fshu 


cl20 


14 18 04.90 


+52 27 40.14 


0.18 


0.07 


24.30 1 






780+!! ,0 ,o 


15+°?? 
0.15 


0.79 


fsh 


cl21 


14 18 05.30 


+52 25 10.63 


0.80 


1.42 


24.44 1 








77+i-? 3 

— 0.77 


0.82 


fsh 


cl22 


14 18 06.51 


+52 33 58.67 


0.84 


0.72 


23. 19 1 








1 29+f- 3 ^ 

1.29 


1.26 


fs 


cl23 


14 18 07.07 


+52 25 23.41 






<26.00 






750+"' 4 

yj. 1 ^o_ ( j qy 


1 78+;-?j 

1 °— 1.00 


2.49 


fsh 


cl24 


14 18 07.33 


+52 30 30.52 


0.85 


0.51 


22.58 1 


0.990 3 


NL 


850+°'°2 


< 13.12 


0.10 


s 


cl25 


14 18 08.06 


+52 27 50.36 






<26.0() 








< 12.52 


1.09 


f 


cl26 


14 1 8 08 94 


4-52 31 50 84 






<26.0() 








339 29+ 271 - 69 


34.55 


f 


cl27 


14 18 09.12 


+52 28 04.04 






<26.00 








6.71+?" 

— 3.67 


6.98 


fshu 


cl28 


14 18 11.26 


+52 30 11.48 


0.86 


1.12 


23. 79 1 


2.910 4 


UNCL 


3 090 +0,22 
o.u»u_q 22 


15 47+45.56 
ld -^' -15.47 


6.54 


fs 


cl29 


14 18 12.16 


+52 28 00.29 






<26.00 








4.00+ 4 ?^ 


3.47 


fs 


<-1 1A 


\A 18 11 1 Q 
14 15 1 9. 1 y 


-US9 31 1 3 Al 






slf. AO 






u.you_ 07 


U * D5 -0.68 


1 if* 
1 . / 


1SI1 


r-1 11 

CI 9 1 


14 16 1 j.jj 


J_S9 1A 1 A QO 


n aa 


1 A^ 

1.4 J 


If. 81 3 
20.61 






1 9^n^ - ^ 

±.2OU_ Qg 


J - U '-1.07 


1 si 
1 .99 


fc 
IS 


cl32 


1 A 1 O IT ild 

14 18 13.96 


i ci o£ o/i to 

+52 26 24.79 






<26.0() 








qo Q'7~l _ 24.61 
'-15.99 


5.17 


fh 


cl33 


14 18 14.27 


i o o i r\ on 

+52 2o 10.99 


1.02 


0.49 


22. 22 1 


2.818 


BL 


2.OO2_ 2 g 


< 1.65 


8.47 


fs 


cl34 


A A AO ICQ*; 

14 18 15.36 


1 CO TO IT £1 

+52 32 4 7.61 


0.00 


1.29 


24. 28 3 








n C7+1-57 
U.8^_ 087 


0.93 


fs 


cl35 


1 A 1 O 1 £. 1CA 

14 18 16.29 


i co on /in Q A 

+52 29 40.30 


0.18 


< 0.01 


20. 07 1 


1.603 3 


BL 


1. 700I 034 


< 0.27 


34.39 


fshu 


^•1 i£ 

Ll JO 


1 A 1 Q 1 (. 
14 16 lO.JJ 


_i_co oc 9/1 ns 
-r jz zj z*+.uj 


A £1 

U.Ol 


n 81 

U.59 


1A 1 8l 
24. 16 








n 58+ 1 - 45 


1 Zl7 

9.4 / 


isn 


cl37 


14 18 16.43 


+52 33 29.77 


0.12 


0.03 


24. 22 1 








< 1.53 


0.29 


fs 


cl38 


14 18 1 (\ 71 
14 16 io. y J 


_I_S9 93 07 QR 


A 75 


A S7 
U.J / 


91 17 1 
Z9.9 / 








s A ^9 


1 91 
1 .Z9 


fsh 


r \ 1Q 

CI jy 


\A 1 C 1 O CiA 
14 16 16.U1 


_i_S9 39 01 98 


n qq 

u.yy 


1 11 

9.9 1 


24 64 1 






n 39A+ 07 


1 ztA+1- 25 
1 -^-0.99 


A 1 7 
u. 1 / 


fs 


C14U 


1 A 1Q 1 Q Q9 
14 16 ly.yL 


j_ co o i i c on 
T" JZ Zl 1J.5U 






s'lf, AO 






l.Ub2_ 12 


q 7fi +16.59 


9 8S 
2.69 


f 
1 


C141 


1 A 1 8 90 ia 
14 15 L\).J\) 


_i_co qq ci no 
-|-jz J J J 1 .Uo 


1 8S 


1 S1 


91 1 Si 
29. 19 






9 9AA^^' 32 
2.2UU_ 32 


n n9 4-4.52 
U.U2_ 02 


S Q7 

D.y 1 


IS 


r-1zl9 
CI 42 


1/118 91 17 
14 15 Zl.J/ 


_l_CO O^ CC AA 
-f-JZ ZD JJ.4-U 






^If. AA 








n 10 
\ O.IO 


1 81 
1 .69 


1 S 


m ai 

CI 49 


1/118 91 1Q 
14 16 21 . JV 


i 19 C/1 qi 

-|-JZ JZ J4-.J1 


n aa 

U.44 


A 9Q 


91 S7l 
29.9 / 








co+4.28 
o.oo_3 22 


A 89 
U.62 


isn 


/.I AA 
CI 44 


1/1 IB Tl "70 

14 15 zi. /y 


i o on cc qo 
+5z zy 5j.oz 


n on 

u.yu 


O f\f. 


1 O 71 1 

iy. / 1 


n ono3 
U.UUU 


^ 1AK 








f.- 


/. 1 A C 


1/1 1 o 99 no 
14 15 22. Uo 


i co o*i en on 
+5z 26 jU. zU 


0.55 


0. 14 


22. 06 1 






0.740I 004 


n cd"f 1 - 31 


A CA 

0.64 


f.- 
tS 


cl46 


14 18 22.41 


4-59 36 07 45 


1.31 


3.83 


24. 15 1 








57+ 1 - 13 


15.50 




cl47 


14 18 22.84 


+52 27 10.11 


1.54 


0.07 


18.26 1 


0.281 5 


AB 


0.530l°;° 4 


< 0.08 


0.04 


fs 


cl48 


14 18 23.07 


+52 21 14.50 


0.00 


0.73 


23.84 3 






1.050lH 6 7 


10.5017-59 


4.81 


fh 


cl49 


14 18 24.97 


+52 23 30.55 


0.66 


0.31 


22.97 1 






1 040+ ' 16 
l.U4U_ Qg 


3.48i 2 ;17 


17.92 


fshu 


cl50 


14 18 25.52 


+52 23 49.48 


0.10 


< 0.01 


22.51 1 










2.52 


fshu 


cl51 


14 18 26.36 


+52 28 18.80 


0.81 


0.68 


23.00 1 








< 0.79 


7.16 


fsh 


cl52 


14 18 26.44 


+52 32 35.01 


0.71 


0.04 


19.91 1 






o.73o+8:gi 


< 0.27 


0.33 


fs 


cl53 


14 18 26.51 


+52 25 59.70 


0.98 


< 0.01 


24.65 2 








12.00+^° 9 f 


5.89 


fsh 


cl54 


14 18 29.76 


+52 27 09.39 


0.56 


0.15 


22.48 1 






0.7201H1 




1.30 


fh 


cl55 


14 18 30.24 


+52 22 12.14 


0.16 


< 0.01 


20.83 1 








< 0.13 


107.14 


fshu 


cl56 


14 18 32.87 


+52 23 49.48 


0.74 


0.06 


20.15 1 








< 2.33 


29.17 


fsh 


cl57 


14 18 37.96 


+52 20 34.62 






<26.00 








< 1.52 


2.53 


s 


cl58 


14 18 38.18 


+52 23 58.55 


1.99 


0.18 


18.84 1 


1.118 5 


BL 


1.500±g;« 


< 0.21 


17.79 


fsh 



